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Abstract 

We analyse recent results on charged particle pseudo-rapidity densities from 
RHIC in the framework of the Dual String Model, in particular when including 
string fusion. The model, in a simple way, agrees with all the existing data 
and is consistent with the presence of the percolation transition to the Quark- 
Gluon Plasma already at the CERN-SPS. 



Recent results on charged particle pseudo-rapidity densities in central Au+Au colli- 
sions, at y/s = 56 and i/i = 130 AGeV, presented by the PHOBOS Collaboration, at 
RHIC, im, give very interesting information that may help to clarify the way the ex- 
pected Quark-Gluon Plasma (QGP) is approached as the energy increases. Those data 
also allow to select among different models of particle production. As in this experi- 
ment the charged particle densities and the average number of participating nucleons are 
simultaneously measured, that provides additional strong constraints to models. 

As nuclei are made up of nucleons, it is natural to start by building nucleus-nucleus 
collisions as resulting from superposition of nudeon-nudeon collisions, in the way it is 
done in the Glauber model approach and generalisations of it. In one (low energy) limit 
the nucleons are seen as structureless and emit particles only in their first collision: this 
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Figure 1: Pseudo-rapidity density normalised per participant pair as a function of cm. energy. The 
lines give predictions for the wounded nucleon model Eq. (|^) (solid line) , the pure multicollision approach 
Eq. (g) (dotted Hne), and the Dual String Model, without fusion Eq. (0) (dash-dotted Hne) and with 
fusion Eq. (|l^ ) (dashed line). A A points are taken from 0, |2^, |l^, pp and pp from ||5[ ^ 0, § 



is the wounded nucleon model P|. The prediction for particle density, when Na nucleons 
from each one of the nuclei in a AA collision participate, is 
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where dN/ dy is the particle rapidity (or pseudo-rapidity) density (for NaNa and nucleon- 
nucleon collisions). If the nucleon is seen as made up of quarks and gluons, with a growing 
number of participating sea quarks and gluons as the energy increases, one anticipates 
dominance of multi-collision processes and the relation 
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to hold, where u^^^ is the number of nucleon-nucleon collisions when A^^ nucleons partic- 
ipate. Elementary multi-scattering arguments |||] give 
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In Fig. H, together with the PHOBOS data, we have presented the quantity ^ 
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as function of the cm. energy a/s for the bounds (|l]) — solid line — and (|2|) with (|) — dotted 



2 



(a) 



(b) 





Figure 2: Two strings are produced by each collisions of valence partons (a) and sea partons (b). Notice 
how the nucleon is broken up in (a) , so that further interactions are of the type illustrated by (b) . See 

the text for further discussion. 



line. We used for dN/dy\,^,^ the parametrisation 0.957 + 0M58\I^{^/s) + 0.0494 ln^(v/i), 



with ^/s in GeV, which fits data from pp and pp non-single-diffractive collisions for cm. 
energies y/s > 22 GeV. The parametrisation used in |[l|, |^ could not be used here because 
it does not fit NA22 data. 

In the Dual String Model (DSM), i.e., the Dual Parton Model [|] with the inclusion 



of strings ||10|| , the limits referred to above appear in a natural way. The valence quarks 
of the nucleon produce particles, via strings, only once — this is the wounded nucleon 
model case — and production is proportional to the number A'^ of participant nucleons 
(Fig^). As the energy and A^i increase the role of sea quarks and gluons increases, they 
interact and produce, again via strings, particles, and the number of collisions z/ becomes 
the relevant parameter (Fig||b). 

One should notice that the diagram of Fig. ||b may be interpreted as multiple inelastic 
scattering, either internally within a given nucleon-nucleon collision or externally involving 
interactions with different nucleons. On the other hand, this diagram may appear repeated 
several times. 

Following 1^, and taking into account the basic diagrams of Fig.s |a and |2|b, we now 
write an expression for the particle pseudo-rapidity density, 
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Na [2 + {2k - l)a] h + {v^ - NA)2kah, 
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where h is the height of the valence- valence rapidity plateau, a is the relative weight of the 
sea-sea (including gluons) plateau and k is the average number of string pairs per collision. 
The diagrams of Fig.s ||a and |]b correspond to /c = 1. However, as we mentioned above. 
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the diagram of Fig. ||b can be iterated with A; > 1 being, in general, a function of energy. 
The number of nucleon-nucleon collisions is, of course. 



and the number Ng of strings is 



= Na[2 + 2{k - 1)] + (z/^^ - NA)2k = 2ku, 
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The first term on the right-hand side of Eq. (U) is just a sum over nucleon-nucleon 
scattering contributions (including internal parton multiple scattering) and we can thus 
write 
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If external multiple scattering is absent, by putting = Na, one obtains the wounded 
nucleon model limit, Eq. ([l|). If multiple scattering dominates, /c ^ 1, we obtain the limit 
of Eq. (1). 

In order to make more transparent the comparison with data, we shall rewrite Eq. (|^, 
by using Eq. (||) and Eq. (|), in the form 



1 dN 



_ dN 

nana 



N 



1/3 



{N]^^ - 1)2(1 - a)h. 
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We show the result of this model in Fig. |T| (dash-dotted line). From comparison of Eq. (||) 
with pp data at low energy, ~ 1, one obtains h ~ 0.75. The parameter a in Eq. (||) was 
put equal to 0.05. 

In the Dual String Model the strings interact, the simplest interaction being fusion due 
to overlap in the transverse plane ||T^. This is the mechanism that leads to percolation 
and to the Quark-Gluon Plasma formation ||Tl|, [l^, |^]. When strings fuse, the strength 
of the colour field is reduced in comparison with the colour field generated by the same 
number of independent strings. This is essentially due to the random sum of colour vectors 
Hi- Qn — ^r=i Qn — V^Q if cill the n strings are of the same type. 

Introducing the dimensionless transverse density percolation parameter rj, 
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where is the string transverse radius (we shall take = 0.2 fm, see ||TT|, [T^), Rna 

1 /3 

the radius of the interaction area {Rna — 1-14A^a ) and A^^ the number of strings, the 
effective reduction factor in particle production is ||T6| , 
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As 7] ^ 0, F{r]) 1 (no fusion) and as oo, F{r]) l/y/v ^ l/^/Wl (all the strings 
fuse) . 

We can consider the parameter rj in two situations. In nucleon-nucleon internal inter- 
actions, and then 
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At present energies rjpp is negligible, r]pp 
external interactions, with 
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10 ^. But we can also consider i] in 
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For A^A 
Eq. 



10^, as in ||l]|, Vn^n^ > ^^Vpp and we shall then only consider rj 



with string fusion becomes 
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In Fig. |l] we have also shown the prediction of the DSM with string fusion (dashed line) 
again with h = 0.75 and a = 0.05. The deviation from the wounded nucleon model limit 
becomes weaker and the agreement with PHOBOS data is quite satisfactory. 
We would like now to make a few comments: 

1. The predictions for particle densities in central Pb+Pb collisions of the DSM 
without fusion and of the DSM with fusion are very different at -y/i = 200 AGeV (RHIC) 
and at y/s = 5.5 ATeV (LHC) as can be seen in the following table, showing the average 
pseudo-rapidity density in the interval [—1,1]: 



cm. energy 


200 AGeV 


5.5 ATeV 


without fusion 


1500 


4400 


with fusion 


700 


1400 



2. The models considered here are essentially soft models. The parameters of the 
elementary collision densities, h and a, were assumed constant, all the energy dependence 
being attributed to the parameter k, the average number of string pairs per elementary 
collision. If h and a are allowed to grow with energy, as a result, for instance, of semi-hard 
effects, the parameter k may then have a slower increase than the one obtained here. 

3. The value found for a, a ~ 0.05, means that the height of the plateau is 
much smaller than the height of the valence- valence plateau. By noticing that for valence- 
valence collisions the two strings stretch all over forward/backward rapidity without much 
overlap, while for sea-sea collisions the two strings do overlap, the value found for a means 



dN 
dy 



0.1 



dN 
dy 



(15) 



val-val 



5 



4. In our Dual String Model with fusion, the parameter Vn^n^ the CERN-SPS has 
the value Vn^n^ ~ 1-8, larger than the critical density {rjc ~ 1.12 -h 1.17) which means 
that percolation transition is already taking place at y/s = 20 AGeV, even allowing for 
non-uniform matter distribution in the nucleus {rjc ~ 1.5) ||r^; this result is valid even 
with k = 1. The observed anomalous J/ip suppression [0 may then be a signature of 



the percolation transition to the Quark-Gluon Plasma 



After the submission of this paper, we became aware of two papers on the same subject 
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1. From the paper of Wang and Gyulassy [T^ we realised that the HIJING point at 
200 AGeV in Ref. |[l|] was 20% too high. This point was corrected in our figure. 

2. In the WA98 Collaboration paper on scaling of particle and transverse energy 
production H^, results on dN/dy were presented in Pb+Pb collisions at 158 A GeV. This 
point was included in Fig. |l| but not taken into account in the calculations. It somewhat 
disagrees with the NA49 point presented in |[l|. 
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